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Abstract. This paper considers the optimization problem of minimizing a 
rational function. We reformulate this problem as polynomial optimization 
by the technique of homogenization. These two problems are shown to be 
equivalent under some generic conditions. The exact Jacobian SDP relaxation 
method proposed by Nie is used to solve the resulting polynomial optimization. 
We also prove that the assumption of nonsingularity in Nie's method can be 
weakened as the finitcness of singularities. Some numerical examples are given 
to illustrate the efficiency of our method. 



1. Introduction 
Consider the problem of minimizing a rational function 

mm 



(1) 



x6i" q(x) 
s.t. hi(x) = ■ ■ ■ = h mx (x) = 0, 
9i(x) > 0, . . .,SVn 2 (» > 0. 



where p(x),q(x),hi(x), gj(x) <E M.[x] := R[xj., . . . , x n ]. As a special case, when 
deg (q) — 0, ([T]) becomes a multivariate polynomial optimization which is NP-hard 
even when p(x) is a nonconvex quadratic polynomial and /ii(x)'s, (?j(a;)'s are linear 

Some approaches using sum-of-squares relaxation to solve the minimization of 
([T]) are proposed in jTOJ [TS] and the core idea therein is in the following. Let S 
be the feasible set of (TTJ). Suppose that r* > — oo and q(x) is nonnegative on S 
(otherwise replace by P ^^" ) ), then 7 € R is a lower bound of r* if and only if 
p(x) — jq(x) > on S. Thus the problem (TIJ can be reformulated as maximizing 7 
such that p(x) — jq(x) is nonnegative on 5 1 , which is related to the representation of 
a nonnegative polynomial on a feasible set defined by several polynomial equalities 
and inequalities. As is well-known, a univariate polynomial is nonnegative on R if 
and only if it can be represented as a sum-of-squares of polynomials (SOS) [25 which 
can be efficiently determined by solving a semidefinite program (SDP) [22l [23] . 
However, when n > 1, due to the fact that a nonnegative multivariate polynomial 
might not be an SOS [25] , the problem (TTJ) becomes very hard even if there are no 
constraints. Denote by M(S) the set of polynomials which can be written as 



mi 7712 



, fi{x)hi{x) + > a j (x)g j (x) 



«=1 .7=0 



where ipi(x) £ WL[x], ga(x) = 1 and <jj(x)'s are SOS. A1(S) is called the quadratic 
module generated by the defining polynomials of S. If M (S) is archimedean, which 
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means that there exists R > such that R— G M(S), then Putinar's Posi- 

tivstellensatz [24] states that if a polynomial f(x) is positive on S, then it belongs 
to M(S). If S is compact, then Schmiidgen's Positivstellensatz [26 states that if a 
polynomial f(x) is positive on S, then it can be represented as 

mi 

(2) f(x) = ^2ipi(x)hi(x) + ^2 av{x)g v (x), 

i=i ue{o,i} m2 

where <fi{x) G M[x], g v {x) = g" 1 ■ ■ ■ g%P 2 2 and a u (x)'s are SOS. The set of polyno- 
mials which have representation (J2|) is called preordering which we denote by P{S). 
Hence, if S in ([1]) is archimedean or compact, we can apply Putinar's Positivstel- 
lensatz or Schmiidgen's Positivstellensatz to maximize 7 such that p(x) — r yq{x) 
belongs to M{S) or P(S). 

In this paper, we present a different way to obtain the minimum r*. Given a 
polynomial / G R[x], let x — (x , x\, . . . , x n ) G ]R n+1 and f hom be the homogeniza- 

tion of /, i.e. f hom (x) — x dcs ^\f{x/xo). We reformulate the minimization of (p} 
by the technique of homogenization as the following polynomial optimization 



(3) 



min p(x) 

s.t. h^ix) = ■■■ = h h ™{x) = 0, q(x) = 1, 
g>?° m (x) >0,...,^° 2 m (i) >0, a; >0, 



where p(x) := x™ x{dcs(p)Acs{q)} p(x/x ) and q(x) := x^ {des(p) ' desiq)} q(x/x ). 
We show that these two problems are equivalent under some generic conditions. As 
a special case, they are always equivalent if there are no constraints in ((T|). The 
relations between the achievabilities of r* and s* are discussed. 

Therefore, the problem of solving ([T]) becomes to efficiently solving problem ((3]). 
For general polynomial optimization, there has been much work on computing the 
infimum of the objective via SOS relaxations, see the survey [TJ] by Laurent and 
the references therein. A standard approach for solving polynomial optimizations 
is the hierarchy of semidefinite programming relaxations proposed by Lasserre [13) . 
Recently, under the assumption that the optimum is achievable, some gradient 
type SOS relaxations are presented in [HI H2] - When the optimum is an asymptotic 
value, we refer to the approaches proposed in [5J H3 [571 HH US] • However, to the best 
knowledge of the authors, the finite convergence of the above methods is unknown 
which means that we need to solve a big number of SDPs until the convergence is 
met. More recently, Nie p~5] proposed a new type SDP relaxation using the minors 
of the Jacobian of the objective and constraints. It is shown [15] that the Jacobian 
SDP relaxation is exact under some generic assumptions. Therefore, in this paper 
we employ the Jacobian SDP relaxation to solve ([3]). 

For another contribution of this paper, we prove that the assumptions under 
which the Jacobian SDP relaxation [T5] is exact can be weakened. Let J be the set 
containing polynomials in the equality constraints and an arbitrary subset of the 
inequality constraints. In order to prove the finite convergence of the Jacobian SDP 
relaxation, it is assumed in [15) that the Jacobian of polynomials in J has full rank at 
any point in V(J) which is the variety defined by J. In other words, if the ideal (J) 
generated by polynomials in J is radical and its codimension equals the number of 
these polynomials, the variety V( J) needs to be nonsingular to guarantee the finite 
convergence. In this paper, we prove that the nonsingularity in the assumptions 
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can be replaced by the finiteness of singularities. More specifically, we show that if 
there are only finite points in V(J) such that the Jacobian of polynomials in J is a 
rank deficient matrix, then the Jacobian SDP relaxation [TS] is still exact. We also 
give an example to illustrate the correctness of our result. 
Another possible and natural reformulation of (TTJ) is 



Clearly, if r* is achievable in (p}, then (j4} is equivalent to ([I]) and we always have 
r* = s*. One might ask why we do not solve ((4]) instead of ([3]). The reason is 
that when we employ Jacobian SDP relaxation [TS] to solve §1 or J|}, we need 
to assume that the optimum is achievable. Actually, s* in (J3j> is more likely to 
be achievable than s* in JJ. To see this, note that when r* is not achievable, s* 
can not be reached either. However, s* might still be achievable when r* is not. 
Some sufficient conditions are given in Theorem 12 . 71 and they are not necessary (see 
Example 14.61 and I4.9|) . For a simple example, consider the problem 



Obviously, r* = s* — and they are not achievable. However, we can reformulate 
it as 



Then s* = and we have two minimizers (0, ±1) which verify that r* is not achiev- 
able by (jcj in Theorem 12.71 

This paper is organized as follows. In Section [21 we reformulate ([TJ) as ([3]) by 
the technique of homogenization and investigate the relations between the achicv- 
abilities of the optima of these two optimizations. We introduce the Jacobian SDP 
relaxation [151 in Section [3] and show that the assumptions therein under which the 
Jacobian SDP relaxation is exact can be weakened. In Section 01 we first return to 
solving the problem (fT]) and make some discussions, then we give some numerical 
examples to illustrate the efficiency of our method. 

Notation. The symbol N (rcsp., R, C) denotes the set of nonnegative integers 
(resp., real numbers, complex numbers). For any t £ 1, |Y| denotes the smallest 
integer not smaller than t. For integer n > 0, [n] denotes the set {1, • • • , n} and for 
a subset J of [n], \J\ denotes its cardinality. For x € R", X{ denotes the i-th com- 
ponent of x. The symbol R[x] = M[xi, . . . , x n ] (resp., C[x] = C[xi, . . . , x n }) denotes 
the ring of polynomials in (xi, . . . , x„) with real (resp. complex) coefficients. For 

a e N", denote |a| = aH \-a n . For x e R™ and a G N", x a denotes x" 1 • • • x£" . 

For a symmetric matrix X, X y (resp., X >- 0) means X is positive semidefinitc 
(resp., positive definite). For u € R™, ||u||2 denotes the standard Euclidean norm. 
C k denotes the class of functions whose k-th derivatives are continuous. 




1 

min — 

xi SK x\ + 
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2. Minimizing Rational Functions by Homogenization 

In this section, we first reformulate the minimization of ([1]) as polynomial op- 
timization by the technique of homogenization and investigate the relations 
between the achievabilities of the optima of these two problems. Then we show 
that the condition under which the problems ([I]) and ([3]) are equivalent is generic. 

2.1. Reformulating the minimization of rational functions by homoge- 
nization. Given a polynomial / £ R[x], let x = (xq,x) — (xq,Xi, . . . ,x n ) £ l n+1 
and f hom (x) be the homogenization of /, i.e. f hom {x) — Xq CS ^ ^ f(x/x ). We define 
the following sets: 

S:={x£ W | h t (x) = 0, gj (x) > 0, i £ [mj, j £ [m 2 ]} , 
(5) So := {Z £ R" +1 | ^ om (£) = 0, gf m (x) > 0, x a > 0, * e [mi], j £ [ma]} , 
S:={x£ R" +1 | h% om {x) = 0, gf m {x) > 0, x Q > 0, * £ [mi], j £ [ma]} • 

Recall that for integer n > 0, [n] denotes the set {1, • • • , n}. Let closure(So) be the 
closure of So in R™ +1 . From the above definition, we immediately have 

Proposition 2.1. f(x) > on S if and only if f hom (x) > on closure(So) ■ 

Proof. We first prove the "if" part. Suppose f m (x) > on closure (So). If there 
exists a point u £ S such that f(u) < 0, then (l,u) G 5q. Thus f hom (l,u) = 
f(u) < which is a contradiction. 

Next we prove the "only if" part. Suppose f(x) > on S and consider a 
point (wo, w ) £ R™ +1 in the closure(S'o). There exists a sequence {(«fc,o,«t)} £ So 
such that lim (uk,o,Uk) — (uq,u). Since Uk t o > for all k £ N, we consider the 

sequence {u/c/ufc,o}- For i = 1, . . . , mi and j = 1, . . . , ma, we have hi(uk/uk.o) = 
^ om (tifc 1 o,«fc)/(«fc,o) deg(hi) = and ffjK/ufc.o) = 5 f m Ko> ^)/K, ) dogfe) > 0. 
It implies that {w/c/ufe,o} €E S*. Thus 



/' lom (u ,u) - lim / hom (u fcl o,«fc) = lim u^ (/) /(« fc /u fc>0 ) > 0, 

fe— >oo fc— >oo ' 

which concludes the proof. □ 



Let d = maxjdeg (p), deg (q)}, p(x) = Xqp(x/xo) and q(x) = x^qix/xo). We 
reformulate the minimization of ([T]) as the following constrained polynomial opti- 
mization: 

r* := min p~(x) 

s.t. h'l om (x) = ■■■ = h%™(x) = 0, q(x) = 1, 
g1 om (x)>0 : ...,gtT(i)>0, x > 0. 

We now investigate the relations between r* and s*. In the following of this paper, 
without loss of generality, we always assume that 

q(x) > on a neighbourhood of a minimizer of ([T]); or q(x) > for 

(6) 

any x £ S with sufficient large Euclidean norm if r is not achievable. 

Otherwise we can replace by ■ Note that we do not assume g(x) is 

nonnegative on the whole feasible set S as in fTU] [JS] . 
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Definition 2.2. [15] If there exists a point ^ (0,u) £ S but (Q,u) £ closure(So), 
then we say S is not closed at oo; otherwise, we say S is closed at oo. 

Theorem 2.3. It always holds that s* < r* , and the equality holds if one of the 
following conditions is satisfied: 

(a) S is closed at oo; 

(b) deg(» > deg((j); 

(c) s* is achievable and Xq > for at least one of its minimizers x* — (xq,x*). 

Proof. We first show that s* < r* . For any u £ S in a neighborhood of a minimizer 

of (TTJ) or with sufficient large Euclidean norm if r* is not achievable, if ^£4 is defined 

at u, then q(u) > by the assumption in ©. Let t — q^) 1 ^ — q{\,u) 1 / d . We 

have q(l/t,u/t) — 1 and (l/t,u/t) £ S, so 

p(l,u) p(l/t,u/t) „ 

. , = —- — - = — — ; rr = P(i/t,u/t) > s , 

q(u) q(l,u) q(l/t,u/t) ^ ' ' 1 1 - 

then we have s* < r* . Therefore, to show r* = s*, we only need to show r* < s* . 

(a) For any feasible point (uo, u) of (J3]), i.e., (uq, u) £ S and <7(ito, u) = 1, since S 
is closed at oo, there exists a sequence {(itfc,o,Wfc)} in S such that itj^o > for any 
k £ N and lim (ufc,o,Ufc) = (uqju). Due to the continuity of q, lim q(v,k,o, Uk) = 1. 

fc— >oo k— >oo 

Hence, we can always assume that for any fe G N, q(u k- o,u k ) > 0. For each 
A; G N, let = q(uk,o, Ufc) 1 ^ and consider the sequence {(wfc.oAfc, Uk/tk)}- We 
have lim (u kfi /t k ,u k /t k ) = {u a ,u k ) and q{u kfi /t k ,u k /t k ) = 1. For i = l,...,mi, 

j = l,...,m 2 , 

^ H^j9j° m (u k ,o,u k ) = g^ m (u kfi /t kl u k /t k ) = -^ufl^g^Uk/u^), 

which imply {u k ,o/t kl Uk/tk) £ S and Uk/ukfi £ 5 for all fc. Hence 

~, x p(uk,o/tk, Uk/tk) p(u k /u k .o) + 

p{u k ,o t k ,u k t k ) = — - — = — — > r 

q(u kfi /t k ,u k /t k ) q{Uk/u k ,o) 

and p(uq,u) = lim p(u k fi/t k ,u k /t k ) > r* which means r* < s* . 

(b) If deg(p) > deg(q), then iro divides q(x). By q(x) = 1, we have uq > for 
any feasible point (uq,u) of ([3]) and it is easy to sec that u/uq £ S, then 

-/ x p(uo,u) p(l,u/u ) p(u/u ) , 

q(uo,u) q(l,u/uo) q{u/u Q ) 

which means r* < s*. 

(c) Since x* > 0, we have x* /xq £ S and 

Py °' ' q(x* ,x*) q(x*/x* ) ~ 
which implies r* = s*. □ 

The following corollary shows that the minimizations of ([1]) and ([3]) are always 
equivalent when there are no constraints in ([1]). 

Corollary 2.4. If mi = = in (H}, then S = R" is closed at oo and r* = s* . 
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Remark 2.5. If S = M™, we can remove xn > in (|3|). In fact, if there are 
no constraints, according to the proof of Part (a) in Theorem \2.3l we only need 
u k,o 7^ to get the same result. Therefore, the global minimization 

■ pO) 

r := mm — — 

iei™ q(x) 

is equivalent to 



(7) 



min p(x) 
s.t. q(x) = 1. 



(8) 



Remark 2.6. We would like to point out that not every S is closed at oo and s* 
might be strictly smaller than r* in this case. For example, consider the following 
problem: 

xi 

■ := mm 

il,i 2 €K (xi — X 2 ) 2 

s.t. xf(xi — X2) = 1, 
x t - 1 > 0. 

Clearly, we have r* = 1. However, [171 Example 5.2 (i)] shows that the set 

{{x l ,x 2 ) £ M 2 I x\{xi - i 2 ) - 1 = 0} 
is not closed at 00. Actually, 

S := {(x 1 ,x 2 ) £ M 2 I x\{x 1 -x 2 ) - 1 = 0,sci - 1 > 0} 
is not closed at 00, either. To see it, we have 

S := {(x ,xi,x 2 ) & K 3 I x\{x\ — x 2 ) — x$ = 0, x\ — x > 0,x a > 0}. 

Consider the point (0, 0, 1) G S. Suppose that there exists a sequence {(xk,o, ^fc.i, x k,2)} 
in S such that lim (xk,o, x k,i, Xk,2) — (0, 0, 1) and Xk.o > for all k € N. Then for 

k— >oo 

< e < 1/2, £/iere exists N € N smc/i t/iai /or any fc > N , we have 
< x fc ,o < e, < £, \%k,2 - 1| < e. 

Thus 

< 4,o = - arjk, 2 ) < 4,i( £ - 1 + e) < 

which is a contradiction. Therefore, S is not closed at 00 and we have s* = < r* 
if we reformulate |[8]) by homogenization as the following problem 

min XqXi 

S.t. (Xl - X 2 ) 2 - 1 = x\{xi - X 2 ) - Xq = 0, 
^1 — Xq > 0, Xo > 0. 

However, in section WJ^ we will show that the closedness at 00 is a generic condition 
for a given set S . 

Let 

S := {x e K™ I hi(x) = 0, gj(x) > 0, i = 1, . . . , mi, j = 1, . . . , m 2 } 

where hi and jjj denote the homogeneous parts of the highest degree of hi and gj , 
respectively. Denote Pd{x) and q<i{x) the homogeneous parts of degree d of p(x) 
and q(x), respectively. 
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Theorem 2.7. If one of the conditions in Theorem \2.3\ holds, then the following 
properties hold. 

(a) r* is achievable if and only if s* is achievable at a minimizer x* = (xq,x*) with 
x* ^0; 

(b) If neither p{x) and q(x) have real common roots inS, norpd{x) and qd{x) have 
real nonzero common roots in S, then s* is achievable. 

(c) If s* is achievable and Xq — for all minimizers x* = (xq,x*) of ([3j) 7 then 
r* is not achievable. For each minimizer x* = (0,x*) of ((3|) 7 if there exists a 
sequence {xk} = {{x ki o,Xk)} in S such that lim x k = x* and Xf. o > for all 

k G N. then lim ^E^Em) = r * . 

Proof. If one of the conditions in Theorem 12.31 holds, we have r* = s*. 

(a) Let x* be a minimizer of then x* G S and t = q{l,x*) 1 ' d = q(x*f/ d > 
by the assumption in ([5]). It is easy to verify that (1/i, x*/t) G S and q(l/t, x*/t) 

1 . We have p(l/t, x*/t) = r* = s* which means (l/t,x*/t) is a minimizer of ([3]). If 
s* is achieved at J* = 6 5 with Xq > 0, then r* is achieved at x* jx% G 5. 

(b) To the contrary, we assume that s* is not achievable. Then there ex- 
ists a sequence {x k } in S such that lim ||ifc|| 2 = oo, lim p(xk) = s* and for 

all fc € N, ?(5fc) = 1. Consider the bounded sequence {5fc/||ifc|| 2 } Q S. By 
Bolzano- Weierstrass Theorem, there exists a subsequence /ll^-l^} such that 
Um Xk-/\\xk- ||o = y for some nonzero y = (yo,y) € S since S is closed. Let 

Pixkj) = s kj , then Jm^Skj = s*. Since p{x k] ) = {\\xk j \\ 2 ) d p(x kj /\\x kj \\ 2 ) and 

lim || = oo, we have p(y) = lim p(xk, /\\xkA\ ) = 0. Similarly, we can prove 

q(y) = lim q{x k - /\\xk- |L) = 0. Thus p(x) and q(x) have real nonzero common root 

y on unit sphere . We have yo = 0, otherwise y/yo'\s& real common root of p(x) 
and g(ar) in S. Therefore = p(y) = p d (y), = q(y) = = ^ om (y) = k{y) 

and < gj° m (y) = 9j{y), i-e. and qd{x) have real nonzero common root y in 

S" which is a contradiction. 

(c) By (jg), if £q = for all minimizers of ([3]), r* is not achievable. Suppose 
x* = (0,x*) is a minimizer of ([3]) and there exists a sequence {aifc} = {(x k fi, x k )} 
in S such that lim x k = x* and x k .o > 0. Then for each k G N, x k /x k .a G S 1 . Since 
p and g are continuous, lim p(xk Oj^fe) = s * an d lim <l( x k,0i x k) = L Therefore, 

k^QG ' fc— >-00 ' 

p(x k /x k ,o) r PO^.o^k) 
lim — = Um —7 — 1 r = s = r . 

fe->oo q[Xk/xk,o) k ^°° q{xk,o,x k ) 

Here completes the proof. □ 

2.2. On the generality of closedness at infinity. Although we have counter 
example in Remark 12.61 we next show that in general a given set S in ([5]) is in- 
deed closed at 00. Therefore, if the constraints in ([1} are generic, (JXJ) and ([3]) are 
equivalent. 

Let us first review some elementary background about resultants and discrim- 
inants. More details can be found in [21 [T7] . Let fi,...,f n be homogeneous 
polynomials in x — (x%, . . . , x n ). The resultant i?es(/i, . . . , f n ) is a polynomial in 
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the coefficients of /i, . . . , /„ satisfying 

Res(fx, ...,/„)= «■ 30 ^ 7i G C", A(«) = • • ■ = /„(«) = 0. 

Let /i,...,/ m be homogeneous polynomials with m < n. The discriminant for 
/i, . . . ,/ m , denoted by A (/i, . . . ,/ m ), is a polynomial in the coefficients of /i, . . . ,/ m 
such that 

A(/i,...,/ w ) = 
if and only if the polynomial system 

f 1 (x)=--- = f m (x) = 

has a solution / a £ C™ such that the Jacobian matrix of /j. , . . . , / n does not 
have full rank. 

We next show that in general a given set 5 in ([5]) is closed at oo. In the following, 
we suppose S is not closed at oo and fix a nonzero point (0, u) G S but (0, u) ^ 
closure^). Let J(u) := {j G [m 2 ] | 5j lom (0,u) = 0}. Then ^ om (0,w) > for 
all j G [ma]\J(u). We have the cardinality mi + | J(u)\ > 1, otherwise (0, u) is an 
interior point of and (0, u) G closure (So). Let 

V(u) := {x G | h^ om (x) = 0, 5 f m (x) = 0, i G [mi], j G J(u)}. 

For any 5 > 0, let 

B((0,«),*) = {(x 0) x) G | \\(x ,x) - (0,«)|| a < 5}- 

Lemma 2.8. Suppose S is not closed at oo, then there exists 5 > suc/i </ia< /or 
<jiZ (xo, x) G -B((0, u), 5) n V(u), we have xq < 0. 

Proof. Suppose such 5 doesn't exist. Consider a sequence {5k} with 5k > and 
lim o"fc = 0. Then for each fc, there exists a point (v,k,o, u k) G -B((0, u), <5fc) n 

such that Ufc^ > 0. By the continuity, there exists N such that for all k > N, 
gj° m (ukfl,Uk) > for each j G [m 2 ]\J(u) which implies (iik,o,Uk) G S for all 
k> N and (0, it) G closure(S'o). The contradiction follows. □ 

Now let us recall the Implicit Function Theorem. 

Theorem 2.9. [12j Theorem 3.3.1; The Implicit Function Theorem] Let 

$(x) = $(xi, . . . ,x n ) = (<f>i(xi, . . . ,x n ), . . . ,cj) m (xi, . . . ,x n )) 

be a mapping of class C k , k > 1, defined on an open set U C M" and taking values 
in R m . Assume that 1 < m < n. Let x° = (x®, . . . , x„) oe a /lied point of U and 
x® = {x\, . . . , Suppose that the Jacobian determinant 

TTien i/iere exists a neighborhood U of x , and open set W C K n_m containing x a , 
and functions /i, . . . , / m of doss C fe on IF smc/i t/iai 

^{xx, . . . ,x n - m , h{x a ), ■ ■ ■ , fm{x a )) = /or every x a G IF. 

Here, x a — (xi, . . . , x n _ m ). Furthermore, f±, . . . , / m are £/ie unique functions sat- 
isfying 

(xeU \ $(x) = j = jx G f7 | x Q G W, x n _ m+fe = fk(x a ) fork = l,...,m\. 
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Let J(u) = {ji 

dh'l 



A(u) 



■iji} and 

— (o,«) ••• 



dh^ 



-(0,«) 



^(0,») 







rai v ' 








SIM • 





Recall that Tij and <?j denote the homogeneous parts of the highest degree of hi 
and respectively. Combining Lemma 12.81 and the Implicit Function Theorem, 
we have 

Lemma 2.10. Suppose S is not closed at oo and mi + \ J(u)\ < n + 1, then 
rank A(u) < mi + | J(u)\. 

Proof. Let m = mi + \J(u)\. Suppose rank A(u) = m. Then there exist m indepen- 
dent columns in A(u). Without loss of generality, we assume the last m columns 
of A(u) are independent, i.e. the Jacobian determinant 



d{ti{ 



horn 



uhom n hom 



n hom\ 

' yji f 



(0,u)^0. 



d\Xfi— 771-1-1 5 • • • , X n *) 

Partition u = (0, u) as (u a , u b ) where u a = (0, u%, . . . , w n _ m ), u b = (u n - m +i, . . . , u n ). 
Then by the Implicit Function Theorem l2.91 there exists an open set W C j^"-" t + 1 
containing u a , and functions fi, . . . , f m of class C k on W such that 

h l l om (x ,...,x n . m J 1 (x a ),...J m (x a )) = 0, i = l,...,m 1 , 

j^K...,^,/!^),-,/™^)^, ie J(u), 

for every z a g W. Here, 5 a = (x , . . - , x n - m ). Therefore, (x a , fi(x a ), f m (x a )) £ 
V(u) for every x a € W. Since W is open and fx, ... , f m are continuous, we can 
choose x a very close to u a such that (x a , fi(x a ), . . . , f m {x a )) £ -B((0, u), S) n V(u) 
with > for any S > 0, which contradicts the conclusion in Lemma 12.81 □ 



The following theorem shows that if the defining polynomials of S are generic, 
then 5 is closed at oo. 



Theorem 2.11. Suppose S is not closed at oo, then 

(a) if mi + |J(u)| > n + 1, then Res{h\ om , h^,gf° m , . 
every . . . , j n - mi +i} Q J{u); 

(b) if mi + \J(u)\ = n, then Res(hi, . . . 



r hom 



) = fo 



(c) if mi + \J(u)\ < n, then A(/ii, 



'mi ) 9ji i 

n ' 9ji j • 



>9j, 



0; 



>9ji 



0. 



Proof. Since h^ om (0,u) = h t (u) = 0, g$ om (0,u) = 9j (u) = for all i € [m x ],j € 
J(m), then the conclusions in (a) and (b) are implied by the proposition of resultants. 
If J(u)\ < n, then by Lemma [2.10[ the Jacobian matrix of (hi, . . . , h mi , gj 1 , . . . 
, gj t ) does not have full rank at u. Hence, the conclusion in (c) follows by the 
proposition of discriminants. □ 
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In this section, we reformulate the minimization of (JTJ as the polynomial opti- 
mization ([3]) by homogenization. Suppose S is closed at oo which is generic and 
always true when S — K", then r* = s*. The relations between the achievabilities 
of r* and s* are discussed in Proposition 12.71 Now the problem becomes how to 
efficiently solve polynomial optimization (J3)). Recently, there has been much work 
on solving polynomial optimization with or without constraints via SOS relaxation. 
In next section, we introduce the Jacobian SDP relaxation [15] and show that the 
assumptions under which the Jacobian SDP relaxation is exact can be weakened. 



3. Jacobian SDP Relaxation Applicable to Finite Real Singularities 
Consider the following polynomial optimization problem 
' fmin ■= min /(x) 

(9) < s.t. hi(x) =--- = h mi {x) = 0, 

gi(x) > 0,...,g m2 (x) > 0. 

where f(x), hi(x), gj{x) € R[xi, . . . , x n ]. In this section, we first introduce the exact 
Jacobian SDP relaxation proposed in [T5] . Then we present our contribution in this 
section by giving a weakened assumption under which the relaxation in [15 is still 
exact. 

Let to = min{mi+TO2, n— 1}. For convenience, denote h(x) — {h\(x) 1 . . . , h mi (x)) 
and g(x) = (gi(x), . . . , g m2 (x)). For a subset J = {ji,...,jk} Q [ma], denote 
gj(x) = (gj 1 (x), . . . , gj k (x)). Symbols Vh(x) and Vgj(x) represent the gradient 
vectors of the polynomials in h(x) and gj(x), respectively. Denote the determinan- 
tal variety of (/, h,gjYs Jacobian being singular by 

Gj = {x e C" | rank B J {x) < mi + | J|} , B J (x) = [V/(x) Vh{x) Vffj(x)] . 

For every J = {ji,..-,jk} C [m 2 ] with k < m~ mi, let rj{, . . . , r}i en {J) bc tlic 
set of defining polynomials for Gj where len(J) is the number of these polynomi- 
als. See [151 Section 2.1] about minimizing the number of defining equations for 
determinantal varieties. For each i = 1, . . . , len(J), define 

(10) ^; 7 (x) = r?/ J] where J c = [m 2 ]\J. 

For simplicity, we list all possible tpf in (|10[) sequentially as 

(11) epi, if2, ip r , where r = len(J). 

JC[m2],| J\<m— mi 

Define the variety 

(12) W := {x e C" | hi(x) = ■■■ = h rni (x) = tpi(x) = ■■■ = ip r (x) = 0}. 
We consider the following optimization 

!f* := min fix) 
s.t. hi(x)=0 (i = l,...,mi), <Pj(x)=0 (.? = l,...,r), 
ff „(x) >0,V^€ {0,1}" 12 , 

where 5^ = .g^ 1 • • • g„T 2 2 . 



MINIMIZING RATIONAL FUNCTIONS BY EXACT JACOBIAN SDP RELAXATION 11 



We now construct N-th order SDP relaxation [13] for ([13]) and its dual problem. 
Let ip(x) be a polynomial with deg {tp) < 2N and define symmetric matrices 
such that 

i>{x)[x] d [x] T d = E A i N)xa > where d = N ~ Tdeg (VO/21. 

QeN":|a|<2Af 

Then the iV-th order localizing moment matrix of ip is defined as 
(14) Lf\y)= £ 

aeN™:|a|<2JV 

where y is a moment vector indexed by a € N n with |a| < 2N. Denote 

L f(y) = E f aVa for ^ = E ^ q:z;Q - 

a6N»:|a|<dcg (/) aeN":|a|<deg (/) 

The AT-th order SDP relaxation [13] for ([13]) is the SDP 

{In : = min L /(y) 
s.t. L < g\ v ) = 0(i = l,...,m 1 ), LW(y) = 0(j = l,...,r), 
LW^O,V^e{0,ir 2 ,yo = l. 

Now we present the dual of ([15]) . Define the truncated preordering generated 
by g-j as 

g(c7„S„) <2Jv' 



piN) = { E ^ 

i/£{0,l} m 2 



ov's are SOS 



and the truncated ideal generated by /i^ and ^ as 



f=l j=l 



deg (r/'ifti) < 27V V« 
deg(<^) < 2iVVj 



It is shown [T3] that the dual of ([T5]) is the following SOS relaxation for (TT 
(16) 



j-(2) 

:= max 7 



i s.t. /(x)- 7 e/ (Ar) + P (Ar) . 

By weak duality we have fffi < fffl < /*. For any subset J = {ji,.--,jk} C [m 2 ], 
let 

V(Mj) = e C" I /K(a?) = 0, a -(a:) = 0, t = l,...,mi, j G J}. 
We make the following assumption. 

Assumption 3.1. (i) mi < n. (ii) Por any feasible point u, at most n — mi of 
gi(u), . . . ,g m2 {u) vanish, (iii) For every J — {ji, . . . , jk} C \m?\ with k < n — mi, 
Jacobian [V/i Vgj] /ias /uH rank on V(h,gj). 

Under the above assumption, the following main result is shown in |15) . 

Theorem 3.2. [TSJ, Theorem 2.3] Suppose Assumvtion \3.1\ holds. Then f* > —00 
and there exists N* G N such that ffP — f& = f* for all N > N* . Furthermore, 
if the minimum f m i n of ([5]) is achievable, then fjy = fffi = f m in f or a ^ N > N* . 



12 FENG GUO, LI WANG, AND GUANGMING ZHOU 

According to Theorem 13.21 it is possible to solve the polynomial optimization 
([9]) exactly by a single SDP relaxation, which was not known in the prior existing 
literature. It is also shown in [TS] that Assumption 13. II is generically true. It is the 
reason why we use this method to solve ([3]) . We will show later in this paper that 
Assumption 13 . 1 1 is always true for ((3]) when the original feasible set S = R™, i.e. for 
the global minimization of a rational function. In the following of this section, we 
prove that the condition (pU| in Assumption 13.11 can always be weakened such that 
the conclusions in Theorem 13.21 still hold. 



Definition 3.3. For every set J — {j±, . . . , jk} C [m2\ with k < n — mi, let 

9j = {x e V{h,gj) | rank [Vh Vg.j] < mi+|J|} and 9 = (J Qj. 

J C [7772],! J\<ri— mi 

We next show that the Jacobian SDP relaxation [TS] is still exact under the 
following weakened assumption: 

Assumption 3.4. (i) mi < n. (ii) For any u € S , at most n—m\ of gi(u), . . . , g m2 (u) 
vanish, (iii) The set 9 is finite. 

Let K be the variety defined by the KKT conditions 

mi 7712 
7=1 j=l 

hi(x) = Hjgj(x) = 0,V(i,i) e [mi] x [m 2 ] 



K={ (x, A, /J.) E £n+mi+m 2 



and 

K x = C™ | (a;, A, /i) € K for some A, ^}. 

Under Assumption 13. 1[ 15., Lemma 3.1] states that W — K x . We now improve this 
result as follows. 

Lemma 3.5 (Revised Version of Lemma 3.1 in |15j). Under conditions (0) and ([m]) 
in Assumption \3.4\ W = 9 U K x . 

Proof. The proof of [HI Lemma 3.1] shows that W\Q QK X CW. With a similar 
argument, we prove 6 C W. Recall that B J = [Vf(x) Vh(x) Vgj(x)]. Choose 
an arbitrary u g 9 and let u £ 9/ for some / C [7712] . If I — 0, then [V/i] and 
B J (u) are both singular for any J C [7712], which implies ipi(u) = and u € W. If 
J 7^ 0, write 7 = {ii, . . . , i t }. Let J = {ji, . . . , jk} Q [mj] be an arbitrary index set 
with mi + fc < m. 

Case I % J At least one j G J c belongs to /. By the choice of I and the definition 

of tpi(x), 

Case I C. J Then [V/i Vpj] and [V/(a;) Wh(x) Wgj(x)] are both singular. 
Hence, all polynomials ipf(x)'s vanish at u. 

Combining the above two cases, we have all (ff(x) vanish at u. Thus, u £ W 
which implies W = 9 U K x . □ 

Lemma 3.6. Under conditions (P and (Jn|) in Assumption \3.4\ if the minimum 
fmin o/® is achievable, then f* = f min . 
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Proof. By the construction of (fT3|) . /* > fmin- Suppose fmin — f{x*) where 
x* is a feasible point of ([9]). If x* ^ O, then the linear independence constraint 
qualification (LICQ) is satisfied at x* which implies x* G K x [501 Theorem 12.1]. 
Since W = 9 U K x by Lemma [3751 we have x* £ W which implies /* = f m in- O 

Next we show that the conclusion in [HI Lemma 3.2] still holds under Assump- 
tion [331 

Lemma 3.7 (Revised Version of Lemma 3.2 in [15 ). Suppose Assumption \3.4\ 
holds. Let T = {x G R™ | gj(x) > 0, j = 1, ...,7112}, Then there exist disjoint 
subvarieties Wq , W\ , . . . , W r of W and distinct V\ , . . . , v r eK such that 

W = Wo U Wi U • • • U W r , W n T = 0, Wi n T ^ 0, i = 1, . . . , r, 

and f(x) is constantly equal to Vi on Wi for i = 1, . . . , r. 

Proof. Denote Zar(K x ) the Zariski closure of K x and let VI = W\Zar(K x ). By 
Lemma I3~5l we have Zar(K x ) C W and £1 C ®. With the proof of [HI Lemma 3.2], 
we can conclude that there exist disjoint subvarieties Wo, W\, . . . ,Wt of Zar(K x ) 
and distinct V\ , . . . , Vt € M such that 

Zar(i^ x ) = Wo U Wi U • • • U W t , W n T = 0, n T ^ 0, i = 1, . . . ,t, 

and /(x) is constantly equal to u, on W, for i = 1, . . . , t. We now consider the set 
Vt. Let Wo = V(E ), then for anyuefin C", T^ U {u} = V(E ) UV((x-u)) = 
V((x — u) ■ Eq). Since ft n C™ C 9 is a finite set by Assumption 13. 4[ if we group 
Wo and SlnC" together then we get a new subvariety. We still denote it by Wo 
for convenience. Then Wo D T = 0. Take any w E VI n E n , if /(w) = w io for 
some io € {1, . . . , <}, then we put w into Wi and get a new subvariety by the same 
reason as Wo- We still write the resulting subvariety as Wi . If for any i 6 {1, . . . , t}, 
/O) ^ u,-, then let W t +i = {w} and v t+1 = f(w) G E. Since fin E™ C 9 is a finite 
set, the above process will terminate and we can obtain the required decomposition 
of W. □ 

Since we get the same result as in [15l Lemma 3.2] under the weakened Assump- 
tion 13. A\ [HI Theorem 3.4] which is based on [HI Lemma 3.2] can be restated as 
follows. 

Theorem 3.8 (Revised Version of Theorem 3.4 in [15]). Suppose Assumption \3.4\ 
holds. Then f* > —00 and there exists N* G N such that for all e > 

(17) f(x)-r+eeI (N " ) +P (N " ) . 

Since e in \17\ is arbitrary, by Lemma 13.61 Theorem 13.21 becomes 

Theorem 3.9 (Revised Version of Theorem 2.3 in [HI)- Suppose Assumption \ 3-4\ 
holds. Then f* > —00 and there exists N* G N such that fffl = fffi = f* 
for all N > N* . Furthermore, if the minimum f m in °f © * s achievable, then 
f { N = f { N = fmin for all N > N* . 

Remark 3.10. We now compare the conditions (iii) in Assumvtion \3.1\ and \3.4\ 
For any J = {ji, ■ ■ ■ , jk} Q H with k < n — mi, suppose the ideal (h,gj) is 
radical and its codimension is m\ + \ J\. Then the condition (iii) in Assumvtion \3. 1\ 
requires the variety V(h,gj) is nonsingular for every subset J. In this section, we 
have proved that if the singularities of V(h, gj) are finite, i.e. the condition (iii) in 
Assumption \3.4\ holds, the Jacobian SDP relaxation [15] is still exact. 
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Corollary 3.11. Suppose that 

(a) For each subset J C [7712] with \ J\ < n — mi, (h,gj) is a radical ideal and its 
codimension is mi + \J\; 

(b) V(h) is a smooth variety of dimension < 2. 

Then the condition (jnT|) in Assumption \3.4\ always holds. Therefore, if conditions 
(0) and (JTTJ) in Assumption \3.4\ are satisfied, then the conclusions of Theorem \3.9\ 
hold. 

Proof. For any subset J C [m 2 ] with \J\ < n — mi, by jlj), 8j is the set of singu- 
larities of V(h,gj). If J = 0, then 6./ = by ©. If J ^ 0, then by 1, Proposition 
3.3.14] or Theorem 5.3], dirnGj < dimV(h,gj). Since dim V(h, gj) < 1 by (jaj) 
and (|b|, Qj is a finite set for each J C [m 2 ] with \J\ < n — mi. Thus the condition 
(jmj in Assumption 13.41 always holds. □ 



We now give an example to illustrate the finite convergence of the Jacobian SDP 
relaxation [15] under the weakened Assumption [ 



Example 3.12. Consider the following polynomial optimization 

min f(xi, x%) := Xix 2 + Xi 

s.t. h(xi,X2) '■= —Xi + x\ = 0. 

Clearly, the minimum f m in — is achieved at (0,0). However, it is easy to verify 
that (0,0) is a singular point and does not satisfy the KKT conditions. Since 
(0, 0) is the only singularity, Assumption 13.41 holds which is also guaranteed by 
Corollary 13.111 In the following, we show the finite convergence of the Jacobian 
SDP relaxation [T5] by giving the exact equation (TIT)) . 

By the construction of (fl~3"|) . m\ — 1,to 2 — and r — 1. ip{xi,x 2 ) '■= 2x2(x\ + 
1) + §x\x 2 . For any e > 0, let 

x 3 

■4>(xi,x 2 ) :=8xi + 8e - Uxfxi - 2Ax\x\ + 24x1x1 + SxixZ + 4x? + 32ex? \ 

e z 

X 1 2x\ X x 2 X\x\ x\ Xi 1 3 2 

8e 3 e 2 4e 3 2 64e 3 8e 2 64e 3 8e 2 
„5^2 j_ 4t ,5 j_ 94^4 _ 243xi°x 2 _ 3x]°X2 _ 45x^a;| ^ Sllxlx 2 
256e 3 !(i 3 128e 3 102 1 3 



+ Ax^Xz + 4x\ + 24xf 



Qt,7„6 o,,3,,4 oq T 9 T 2 90^.6^,2 zI^t- 4 ^ 4 4 ^.r. 10 <-r 4 
0x^2 1 J-t' ^ 2 ^ ojx^2 Zjc?X2o-2 ^uJj^2 ^tuu-^ 0^2 



32e 3 32e 2 8e 2 32e 2 1024e 3 64e 3 

r 6 
' 2 



\*Ix^X<2 3x^3^2 3x^X2 X^X<2 3x^X2 

32e 2 + 2e 2 + 4e 2 + 1024e 3 ~ 256e 3 ' 



^l,a-2j •— orl _ a 100.3 + erf_3 C10.3 ZX 1 X 



™ 10 ™5 
Xi x 2 


15x}°x 


3 /y»9 ™3 

2 1 x l x 2 


7 5 

x l x 2 


81a;}°x 2 


32e 3 


128e 3 


4e 2 


64e 3 


512e 3 


15x-^ X2 


- 4x^x 2 


«,6^,3 
1 x l- i 2 


45 

x l x 2 


337x^X2 


256e 3 


8e 2 


512e 3 


2048e 3 


15x-^ X2 


- 2x^X2 


™3 ~,3 
, • i l x 2 


4 

x 1 a; 2 


lx\x 2 , 


2048e 3 


64e 2 


16e 3 


We 2 


xix 2 
128e 3 


X2 

16e 2 


4ex 2 . 







\\x\x2 

2+ 16e 2 
, 7 59x?x 2 
'+4^2 + ^2- 

2x\x2 — 4xiX2 
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ao(xi,x 2 ) :=16 e + 



(xiXg + X\ + l) 2 

4 





xia; 2 + x\ (x\x\ 



2s 



It can be verified that 

f(xx,x 2 ) + e = o- (x 1 ,x 2 ) + ip(x 1 ,x 2 )h(xi,x 2 ) + 4>{xi,x 2 )ip(x 1 ,x 2 ). 

Since each term on the right side of the above equation has degree < 20, we take 
N* = 10 in (|17[) . Because o~q{xi,x 2 ) is a sum of squares of polynomials, we have 
cr (xi,x 2 ) £ p( 10 ) and 4>(x\, x 2 )h(x\, x 2 ) + 4>{x\, x 2 )ip(x\, x 2 ) € I^ w \ Therefore, 
f(xi,x 2 ) +ee + P( 10 ' for any e > 0. Hence, we have ffl = f^ = f min = 
for all TV > 10. □ 

A practical issue in applications is how to detect whether (|15l) is exact for a 
given N. Nie [T^ pointed out that it would be possible to apply the flat- extension 
condition (FEC) |S]. When FEC holds, (|15l) is exact for © and a very nice software 
GloptiPoly [9 provides routines for finding minimizers if FEC holds. In general, 
the FEC is a sufficient but not necessary condition for checking finite convergence 
of Lasserre's hierarchy. More recently, Nie |16j proposed the flat truncation as a 
general certificate. For the polynomial optimization (|9]), define 



When ip = 1, V\ (y) in (|14p is called moment matrix and is denoted as Mjy(y) := 
L^\y). For a given integer JV € N, we say an optimizer y* of (TT5]) has a flat 



Assuming the set of global minimizers is nonempty and finite, [161 Theorem 2.2 and 
2.6] show that the Putinar type or Schmiidgen type Lasserre's hierarchy has finite 
convergence if and only if the flat truncation holds. As an application, [TBJ Corollary 
4.2] also points out that if ^ has a nonempty set of finitely many global minimizers 
and Assumption 13 . 1 1 is satisfied, then the flat truncation is always satisfied for the 
hierarchy of Jacobian SDP relaxations. Since we have proved that Assumption 13. II 
can be weakened as Assumption 13.41 we have 

Corollary 3.13 (Revised Version of Corollary 4.2 in [H]). Suppose © has a 
nonempty set of finitely many global minimizers and Assumption \3.4\ is satisfied. 
Then, for all N big enough, the optimal value of (|16p equals the global minimum of 
@ and every minimizer of (|15[) has a flat truncation. 



In this section, we return to the minimization of ([TJ . We first apply the Jacobian 
SDP relaxation discussed in Section [3] to reformulate §5§ as (fl"3"| for which we 
consider the finite convergence of the SDP relaxations. Next, we do some numerical 
experiments to show the efficiency of our method. 



d h ,i = \deg(hi)/2] , d gd = \deg(g 
d = max{l, dh,i, ■ • ■ , <4,mi , d g ,i, ■ ■ . 





4. Revisiting Minimization of Rational Functions 
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4.1. Minimizing Rational Functions by Jacobian SDP Relaxation. Con- 
sider the number of new constraints added when we employ Jacobian SDP relax- 
ation to solve ((31) • As mentioned in |15j . the number of new constraints in (|13[) is 
exponential in the number of inequality constraints. Hence, if the number of in- 
equality constraints is large, (|13[) becomes more difficult to solve numerically. In the 
following, we employ the Jacobian SDP relaxation to reformulate ((3|) as (|13| . We 
show that the number of the new equality constraints ip^s in (|13j) can be reduced 
due to the special inequality constraint x$ > in ©. 
In ([3]), for convenience, let 

-=q(x)-l = 0, g^™i(x):=x Q >0 and m := min{mi +m 2 + 2, n}. 

Denote 

V (— — —\ 

\dxQ 1 dx\ ' ' dx n J 

According to (jlOl) and (TTT1) . we need to consider all subsets of [m2 -I- 1] with 
cardinality < m — mi — 1. Let I = min{m — mi ~ 1,7712}. We first consider 
the subsets without mi + 1, i.e., every subset J = {ji, ■ ■ ■ ,jk} Q [ m 2\ with 
k < I. Denote h hom = {h\ om , h>>°™, ft^+i) and 9 J° m = (5ji° m > • • -i9j° m )- Let 
{771, . . . , r)i en (j)} be the set of the defining equations for the determinantal variety 

Gj := {£ € C ,l+1 I rank [Vsp V^' lom V £ ^ om ] < m x + \J\ + 1}. 
For each i = 1, • • • , len(J), define 

^f(i) = % II 9j 0m (x), where J c = [m 2 + 1]\J. 

For every subset J considered above, denote J' = JU {7712 + 1} C [7712 + 1]. It can 
be checked that the collection of these J's and J"s contains all subsets of [7772 + 1] 
with cardinality < m — mi — 1 and some possible J"s with cardinality = m — mi 
(which will happen when 77 < mi + 7772 + 2). 

Case I J' I < 777 — mi — 1 All these J"s compose of the subsets of [7772 + 1] containing 
7772 + 1 with cardinality < m — mi — 1. It is easy to see that the set of the defining 
equations for the determinantal variety 

G r := {Z € C n+1 | rank [V £ p \7 £ h hom V ig h j om V £ x ] < mi + \J\ + 2} 

is a subset of {771, . . . , ?7/ e „(j)}. We generally suppose it to be {771, . . . , 77t(j)} with 
t(J) < len(J). For i = 1, ••• ,t(J), define 

^{(x) = r ]l - Y[g!} om (x), where J' c = [m 2 + 1]\ J 1 . 

jeJ' c 

Case I J'| = m — mi It is easy to check that Gji = C" +1 . Thus for convenience, 
we set t(J) — in this case. 

Then for every subset J C [m 2 ] with | J| < I, we have 



(18) 



<Pi(x) = fi (i) • x , i = !,••• ,t(J). 
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Now consider the SDP relaxations [13] for the following polynomial optimization 



' p* :— min p(x) 

iGR™+ 1 

s.t. h\° m {x) = ■■■ = h%™(x) = h h ™ +l {x) = 0, 



(19) 



< 



<pi{x)=<S,<p J j (5) = 

(i € [len(J)],j E [t(J)],J C [m 2 ], |J| < ?), 

s£ om (z) > 0,V^G {0,l} m2+1 , 



where g£ om = {g^ om ) Vl . . . (g^J^+i. We now show that for each J C [m 2 ] 
with | J | < Z, constraints ip{ (x) = ■ ■ ■ = ipf,j*.(x) — can be removed from (T19]) . 
Consider the N-th order SDP relaxation ([T5|) for ([T9")) . By ([JS]) and the properties 
of localizing moment matrices in |141 Lemma 4.1], we have 



L^)(y) = implies L^'(y)=0, j = 1,. . t(J), J C [m 2 ], \J\ < I. 
In the dual problem (fT5|) . by (IT51) . the truncated ideal 

(len(J) t(J) \ mi+X 1 

E kri +Y,Q<pf + E ^ om f > where 



Vi,j,*, deg < 2JV,deg(0v/') < 2iV, deg (^ fc ^ om ) < 27V, 



agrees with 
(20) 



mi+l 



E 

JC[m 2 ],|J|<Z Vi=*(J)+l 



Zen(J) t(J) 

E hvi +Y,tj<pf I + E ^ om } where 



fc=i 



V»,i,A, deg(<^f) < 2AT,deg (<>/') < 2N,deg(^ k h h k om ) < 2N. 

Therefore, we can remove tp{(x) — ■■■ — = in (I19[) and improve the 

numerical performance in practice. Hence we consider the following optimization 

f p* := min p(x) 

s.t. h h x ° m {x) =■■■ = h h ™(x) = - 0, 

< 2ii { =0,<pf(x) =0 

(i = i(J) + l,...,Zen(J),je [t(J)],JC [m 2 ],|J|<0. 
.9^ om (i)>0,V^e{0,ir 2+1 , 
The JV-th order SDP relaxation [13] for dHJ) is the SDP 
' p%' := min Lp(y) 



(22) 



s.t. L^L{y) = ■ ■ ■ = L^Uy) ^ L^L (y)=0, 



L$(y) = 0,L^(y) = Q 
{i = t(J) + l,...,len(J),je[t(J)],JC[m 2 ],\J\<l), 
L^L ^0,V^G{0,l} m2+1 ,yo = l. 



-W, 
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The dual problem of fl22J) is 

(23) pP ■= max i 7 s.t. p(x) - 7 G 7 (A,) + P w . 



7GR" +1 

where is the ideal defined in (l20l) and 



_;ve{0,l} m 2 + 1 



^ Jiora 



deg (ov<?£ om ) < 2iV 
av's are SOS 



Definition 4.1. For every set J — {ji, . . . , jk} C [777,2 + 1] with k < n — mi, let 
e, 7 = {ie V(h hom ,g h j om ) I ran* [V 5 /i' lom V s c$ oro ] < mi + \J\ + 1} 

and 

e = |J Qj. 

JC[m 2 +l], |J|<ra-mi 

Assumption 4.2. (i) mi < 77; (ii) For any u E S in (JSJ) ^ at most n — mi of 
5 5 tom (u), . . . ,g^i(u) vanish; (hi) T/ie sei 9 zs /sra«e. 

By Theorem 12.31 and 13.91 we have 

Theorem 4.3. Suppose Assumption \4-2\ holds. Then p* > —00 in (|2ip and i/iere 
exists N* G N suca that p^r = pff = p* for all N > N* . Furthermore, if one of 
the conditions in Theorem \2.3\ holds and the minimum s* of ([3J is achievable, then 

p w = p p = r * j or all N > iV *_ 

Corollary 4.4. If S = R" in ([1]) and s* is achievable in (O, t/ien £/iere exists 
N* G N swc/i </iat p$ = = r* /or a// N > N* in ([22} and d23j . 



Proof. Since g is homogeneous, regarding and x as vectors in then d- q = 

Vq T ■ x by Euler's Formula. Thus V(q - 1) = Vq = implies £ = 0, i.e. 6 = 0. 
Hence, Assumption 14.21 is always true for ([7]). Then by Corollary 12.41 and Theorem 
14.31 the conclusion follows. □ 



In the end of this subsection, we would like to point out that s* in ^ might 
not be achievable in some cases. If the inhmum of a constrained polynomial opti- 
mization is asymptotic value, some approaches are proposed in [5l 129). Hence, we 
can still use these approaches to solve ((3]). However, to the best knowledge of the 
authors, the hnite convergence for these methods is unknown. 



4.2. Numerical Experiments. In this subsection, we present some numerical 
examples to illustrate the efficiency of our method for solving minimization of |T]) . 
We use the software GloptiPoly 9 to solve (|22} and ([23]). 

4.2.1. Unconstrained rational optimization. In the following, Example 14.51 and 14.61 
are constructed from the Motzkin polynomial 

(24) M(xi, £2, £3) = xfx2 + x\x\ + x\ — ix\x\x\. 

As is well-known, M(x\, xi, £3) is nonnegative on ]R 3 but not SOS [25] . 



MINIMIZING RATIONAL FUNCTIONS BY EXACT JACOBIAN SDP RELAXATION 19 



Example 4.5. 18, Example 2.9] Consider minimization 
25) mm r(x 1 ,x 2 ) := ? V . 

Taking = 1 in Motzkin polynomial, we have + 1 — 3x^1 ^ on R 2 . 

Since r(l, 1) = 3, we have r* = 3 and there are four minimizers (±1, ±1). However, 
x\x\ + x\x\ + 1 — r*x\x\ is not SOS. To solve this problem, the authors in [18] 
used the generalized big ball technique. More specifically, it is assumed that one of 
the minimizers of (|25[) lies in a ball B(c,p) and the numerator and denominator of 
r(xi, X2) have no common real roots on B(c, p). However, it is not easy in general 
to determine the radius p of this ball. We now solve this problem using our method 
without the assumptions in [18] . 

We first reformulate the problem as the following polynomial optimization: 

min p(xo, X\, X2) ■= x\x\ + x\x\ + .Tq 

Xq,Xi ,X2 

s.t. q(xo, xi, X2) := x\x\x\ = 1. 

By Jacobian SDP relaxation (|2~lj) . we need 3 more equations: 

ipi(x ,xi 1 x 2 ) = Axlxlxl(xl -x%) = 

<P2(xo, x\,X2) = ^x\x\xo{2x\x\ + x\x\ — 3xq) = 

(f3{xo, xi,x 2 ) = Ax\x 2 xq{x\x\ + 2x\x\ — 3xq) = 

By the condition x\x\x\ = 1, the above three equations can be simplified as 

tp 1 (x ,x 1 ,x 2 ) = x\ - x\ = 
(p 2 (xo, xi,x 2 ) = 2x\x\ + x\x\ — 32)g = 
^3(^0, xi, x 2 ) — x\x\ + 2x\x% — ix% = 
We need to solve the following new problem 
min x\x\ + x\x\ + x® 

s.t. x-^X2Xq 1 — 0, c 2^x^X2 ~\~ x^x^ 3^/Q — 0, 

X-^ X*2 — 0, X-yX2 ~\~ 2x^3^2 *^^0 ~~ ^' 

Using GloptiPoly to solve this problem, we get the following results: 

• N = 3. The optimum is 3, but extracting global optimal solutions fails. 

• N = 4. We get 8 optimal solutions for (xq, x\, x 2 )- (±1, ±1, ±1) from which 
we get all the optimal solutions for original problem: (±1, ±1). □ 

Example 4.6. [TSl Example 2.10] Consider the following problem 

p(xi,x 2 ) _ x\ + x\ + x% 



(26) min r(x\ , X2) ■ 



xi,x 2 eM q(x\,X2) 



Taking x 2 — 1 in Motzkin polynomial (|24j) . we have r* — 3 with 4 minimizers 
(±1,±1). The denominator and numerator have real common root (0,0). In |18) . 
the SOS relaxation extracts 6 solutions, 2 of which are not global minimizers but 
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the common roots of p{x) and q(x). We reformulate it as the following polynomial 
optimization and solve it by Jacobian SDP relaxation. 



(27) 



min p(xo,xi,X2) '■= XjXq + x^Xq + x? 

10,11,126! 

2„2 2 



s.t. q(xo,Xi,x 2 ) := x 1 x 2 x = 1. 



Using GloptiPoly, we can still extract 8 solutions of (|2"T|) and obtain all the 4 optimal 
solutions of (|26|) as in Example 14.51 In our method, the constraint q(xo, x\, x 2 ) = 1 
prevents extracting the common real roots of p(x) and q(x). This example also 
shows that Condition (jbj in Theorem 12. 71 is only sufficient but not necessary. □ 

The following example is generated from the Robinson polynomial 

6 , 61 6 _ ( 4 2 , 24, 42, 24, 42, 2 4\ , o 2 2 2 

which is nonnegative on R 3 but not SOS .25]. 
Example 4.7. Consider the following problem 

p(x 1 ,x 2 ) _ x\ + x\ + ?>x\x\ + 1 



(28) min r(xi,x 2 ) .— - — — , , -, — , , , 

x u x 2 eM v ; q(xi,x 2 ) x\{xl + l) + x\{x\ + l) + (x\+xf) 

Taking X3 — 1 in Robinson polynomial, we have p(x\ , x 2 ) — q(x\ ,x 2 ) > on 
Since r(l, 1) = 1, r* = 1. We reformulate it as the following optimization: 



(29) 



min x\ + x 2 + 3x 1 x 2 x + Xq 



S.t. X-^(x 2 -\- Xq) "i" X2(x^ -|- Xq) ± Xq(x^ ± X2) — 1- 

The numerical results we obtained are: 

• For relaxation order N — 5,6, we get the optimum s* = 1, but the mini- 
mizers can not be extracted. 

• For relaxation order N — 7, we extract 20 approximate minimizers of (|29p : 

(-0.0000, ±0.8909, ±0.8909), (±0.8909, ±0.8909, -0.0000), 
(±0.8909, -0.0000, ±0.8909), (±0.7418, ±0.7418, ±0.7418). 
The above solutions correspond to the exact minimizers of (f29|) : 

(0, ± ^,±^), (±£.±£.<». 
1 11 1 1 

(± ?1' ' ± W' ^W^^W' 

There are four solutions with the first coordinate Xq = which indicate that mini- 
mum r* — 1 is also an asymptotic value at 00 by Theorem 12.71 In fact, 

p(xi,x 2 ) 
hm — = 1 = r . 

2:1,3:2^00 q[xi, x 2 ) 

From the other 16 solutions, according to (jaj in Theorem 12. 7\ we get 8 global 
minimizers of ([28]): (±1,±1), (±1,0), (0,±1). □ 

Example 4.8. [18l Example 3.4] Suppose function ip(z) and cj)(z) are monic com- 
plex univariate polynomials of degree m such that: 



ip(z) = z m + i> m -\z m + ■ ■ ■ + ipiz + Vi 
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It is shown in 11 that finding nearest GCDs becomes the following global mini- 
mization of rational functions 

p(xi,x 2 ) _ \ip(xi + ix 2 )\ 2 + \<j>(xi + ix 2 )\ 2 



(30) 

__ __ -/ ~ «... a m — 1 



k=a 

where deg(p) = 2m and deg(q) = 2(m — 1). Let 

V>(z) = z 3 + z 2 - 2, 0(z) = z 3 + l.bz 2 + 1.5z - 1.25. 

Using our method for relaxation order N = 5, we get four optimal solutions of the 
optimization reformulated from (1301) by homogenization: 

(0.7050, -0.7073, ±0.7763), (-0.7050, 0.7073, ±0.7763). 

The corresponding minimizers of (|30[) are 

(x! w -1.0033, x 2 w ±1.1011), (ari w -1.0033, x 2 ~ ±1.1011) 

which are the same as in [18j . The minimum is r* 0.0643. □ 

4.2.2. Constrained rational optimization. We now give some numerical examples of 
minimizing of rational functions with polynomial inequality constraints. We first 
consider an example for which p{x) and q(x) have common roots. 

Example 4.9. [TS] Consider the following optimization 

(31) min r(x) := + * s.t. (1 - x 2 ) 3 > 0. 

V ' x€SL (1 - X 2 ) 2 

As shown in [18] . the global minimum r* — || w 0.8438 and the minimizer 
x* = — -| rs —0.3333. If the denominator and numerator have common roots, SOS 
relaxation method proposed in |18j can not guarantee to converge to the minimum. 
Reformulating the above problem by homogenization, we get 

min Xq + xiXq 

(32) \ s.t.4-2^l4 = l, 

Xq — 3xqX 2 + SxqXi — x\ > 0, xq > 0. 

For relaxation order N — 7, by the Jacobian SDP relaxation, we get the optimal 
solution of x* w (1.0607,-0.3536) and the minimum s* « 0.8437. According 
to in Theorem 12.71 we find the minimizer of (|3"Tj) : ps —0.3334. □ 



We next consider Example 14.71 with some constraints. 
Example 4.10. Consider optimization 

/qq\ * . p(xi,x 2 ) x\ + x\ + 2,x\x\ + 1 



xes q{ Xl ,x 2 ) x\{x\ + 1) + x 2 (xj + 1) + (xf + x\)' 

(a) S — {{x\,x 2 ) e R 2 | x\ + x\ < 1}. It is easy to check that S is closed at oo. 
By Theorem 12.31 r* is equal to the optimum of the following optimization: 



min x®+X2+ 3xfx 2 xl + Xq 

S.t. X-^(x2 ± ^0) ± ^2(^1 ± "^o) ± ^o(^l ± ^2) = 

— x\ — x\ > 0, XQ > 0. 
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For relaxation order N = 7, we get r* = s* = 1 with 4 approximate minimizers: 
(0.8909, ±0.8909, -0.0000), (0.8909, -0.0000, ±0.8909), 

which correspond to the exact minimizers: 

1 11 1 

Then we get four minimizers of ([55)1 : (±1,0), (0, ±1). 
(b) S = B(0, V2) c = {(^l,^) e K 2 | x\ + x\ > 2}. S is noncompact but closed at 
oo. By Theorem 12. 31 we solve the following equivalent optimization: 



s* := min x\ + x\ + Zx\x\x\ + Xq 

S.t. X-y{x<2 ± *^o) ± X<2 (x ] ± *^o) *^o(*^l ^2) — 

xj + xj- 2x1 > °> > 0. 
For relaxation order N — 7 , we get r* = s* = 1 with 8 approximate minimizers: 
(0.0002, ±0.8909, ±0.8909), (0.7418, ±0.7419, ±0.7419), 

which correspond to the exact minimizers: 

1 1 11 1 

( ' W ± W' ( ^' ± ^' ± ^' 

The former solutions indicate that r* = 1 is also an asymptotic values at 00. 
From the latter solutions, we get four minimizers of (|33l) : (±1, ±1). □ 
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